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We introduce a simple lattice model with Ising spins to explain recent experimental results on 
spin freezing in a hollandite-type structure. We argue that geometrical frustration of the lattice 
in combination with nearest-neighbour antiferromagnetic (AFM) interactions is responsible for the 
appearance of a spin-glass phase in presence of disorder. We investigate this system numerically 
using parallel tempering. The model reproduces the magnetic behaviour of oxides with hollandite 
structure, such as a — Mn02 and presents a rich phenomenology: in absence of disorder three 
types of ground states are possible, depending on the relative strength of the interactions, namely 
AFM ordered and two different disordered, macroscopically degenerate families of ground states. 
Remarkably, for sets of AFM couplings having an AFM ground state in the clean system, there 
exists a critical value of the disorder for which the ground state is replaced by a spin-glass phase 
while maintaining all couplings AFM. To the best of our knowledge this is the only existing model 
that presents this kind of transition with short-range AFM interactions. We argue that this model 
could be useful to understand the relation between AFM coupling, disorder and the appearance of 
a spin- glass phase. 
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I. INTRODUCTION 



Spin-glasses are magnetic phases where disorder and 
frustration suppress any simple ordered patterns, like 
ferro- or antiferromagnetic states, and have instead a 
spin- freezing transition into an amorphous glassy ordered 
state at low temperatures (or other control parameters). 
Such glassy states feature many interesting and unusual 
properties, like power-law correlations in the absence 
of any broken symmetry 1 and non-trivial long-time be- 
haviour. The free-energy landscape of such systems is 
very rough and has many metastable states that are sep- 
arated by high barriers 1 . 

The two crucial ingredients necessary to produce a 
spin-glass are disorder and frustration. The canon- 
ical representatives of spin-glasses, magnetic alloys, 
have an oscillating long-range Ruderman-Kittel-Kasuya- 
Yosida (RKKY) spin-spin interaction 5 . The stan- 
dard Edwards- Anderson model has quenched random 
nearest-neighbour couplings of both signs which mimics 
the frustrating effect of the longer-ranged RKKY inter- 
action. In both cases, as well as in many others, it is 
the interaction, which is both disordered and frustrating 
{i.e. mixture of ferromagnetic (FM) and antiferromag- 
netic (AFM) couplings) that is responsible for appear- 
ance of a spin-glass. There is also a large class of materi- 
als where frustration has geometrical origin: the combi- 
nation of the AFM interactions and geometry of the lat- 
tice suppresses the natural AFM order 5 and makes the 
system extremely susceptible to perturbations. In this 
case often even a small disorder in the coupling strengths, 
that does not change their AFM character, is enough to 
obtain a spin-glass (> ~'\ Consider now a geometrical frus- 



trated system with nearest-neighbour AFM couplings of 
different strengths. Can the AFM ground state still sur- 
vive? If so, is it possible to obtain a spin-glass phase 
by tuning the disorder and maintaining the AFM inter- 
actions? Do we need to introduce a small amount or a 
large disorder to break the AFM ground state? 

In this paper we introduce a model which addresses 
these points and extends the previously studied models 
of spin-glass transition in presence of geometrical frus- 
tration. The interactions in the model arc geometrically 
frustrated with all the couplings being AFM but of differ- 
ent strengths. The natural AFM order is not always sup- 
pressed in the clean case or by disorder; a critical amount 
of disorder is required to suppress it. We can tune the 
strength of the geometrical frustration by changing the 
relative values of the couplings and obtain three different 
types of ground states: one is non-degenerate and AFM 
ordered; the other two arc exponentially degenerate and 
disordered. The disorder in the coupling strengths which 
leaves the couplings AFM, puts all these ground states in 
competition and leads to the spin-glass state. This pro- 
vides an example of a system where the geometry forces 
different ground states to compete and there is a direct 
transition from an AFM state to a spin-glass phase by 
tuning the disorder, while having only AFM couplings. 
To the best of our knowledge this is the only existing 
model that presents this kind of transition with short- 
range, nearest-neighbour AFM interactions. 

This model has been developed to explain the magnetic 
behaviour of manganese (Mn) oxides with the hollandite 
structure, such as a-Mn0 2 "' 1 ", Ki. 5 (H30) K Mn 8 0i 6 u , 
Bai.2MngOi6 12 , etc.. As we will show below, experimen- 
tally these materials have a rich magnetic phenomenology 
whose origin is poorly understood. In particular it is not 



clear if the magnetic properties can be explained by con- 
sidering only the bulk properties like the effects of doping 
and frustration or if surface effects have to be taken into 
account. Therefore a theoretical work is needed for a 
better understanding of the origin of the magnetic be- 
haviour of these materials. This kind of studies will open 
the venue for a systematic approach to tune and optimise 
the properties of these compounds. 

Hollandite-type Mn oxides are nanoporous materi- 
als 9 ' 10 that are interesting on its own for their large 
number of applications, as ionic conductor, as catalysts 
for oxygen reduction and oxygen evolution reactions, in 
lithium-air batteries 1 ' 5-15 , in supercapacitors 16 , for the 
energy extraction from salinity differences 1 ' and as a 
water oxidation catalyst 18 . The large lateral size of 
the channels (see Fig. 1 (a)), of the order of 0.46 nm, 
make it possible for some big cations such as K + , Na + , 
Ba 2+ and H3 + to be introduced during material syn- 
thesis '', thus opening the possibility to tune the 
magnetic, chemical and physical properties by cation 
doping. The magnetic nature of this compound origi- 
nates from the presence of the manganese ions, which 
have localised magnetic moments due to the open shell 
of d electronic states 11,10 . These properties together with 
the presence of magnetic frustration present because of 
the existence of structural triangles (see Fig. 1 (c)), are at 
the origin of the rich variety of magnetic ordered phases 
that have been identified experimentally in these materi- 
als. 

For example a-Mn02 in absence of doping elements in 
the channels, i.e. in absence of disorder, has an AFM 
phase at low temperatures ' . Luo et al. have investi- 
gated the dependence of the magnetic phase on the dop- 
ing of potassium ions 20 > 21 ) and have shown that above a 
critical concentration of K + ions, a difference in the mag- 
netic susceptibilities between zero-field-cooled and field- 
cooled samples is observed which is interpreted as the 
onset of the spin-glass phase. When all the sites avail- 
able for K insertion became full, the AFM behaviour 
was recovered. The spin-glass phase was previously ob- 
served in the case of the KMngOi6 material 1 " with the 
same behaviour of the magnetic susceptibilities. Lan and 
co-workers proposed a surface effect to be responsible 
for the spin-glass behaviour in sodium-doped a-MnC>2 
nanorods 2 ' 5 . They proposed that, at low temperature, 
the core of the rod is in an AFM state while the sur- 
face spins contribute to a net magnetic moment of the 
specimen 21 . Their explanation thus relies on the high 
surface area present in the channels of the material more 
than on the doping. The magnetic behaviour can be- 
come even more complex by including more than one 
doping species. For example Sato et al. have inves- 
tigated Ki. 5 (H 3 O) a; Mn 8 Oi 6 (0<x<0.5), and found three 
magnetic transitions by varying the temperature 1 1 . They 
explained their results on the basis of a FM phase exist- 
ing at intermediate temperature, and a ground state with 
helical magnetism appearing below 20 K. 

Due to the complex magnetism present in these mate- 
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Figure 1. (Colour online) Different views of the hollan- 
dite structure for the case of a-Mn02 compound 20 ' 21 (a) 
the xz plane containing both Mn and oxygen (O) atoms (b) 
panoramic view of the channel with only Mn atoms (c) Struc- 
tural triangles composed by magnetic atoms that play an im- 
portant role in setting the magnetic properties of the oxide. 



rials, it is important, as a first step, to disentangle the 
contributions of bulk effects (doping and frustration) and 
surface effects in order to gain a better understanding and 
control of the properties of these compounds. The model 
introduced in this paper is able to reproduce the essential 
features of these systems, in particular the transition tak- 
ing place at increasing doping/disorder between an AFM 
ground state and a spin-glass phase, on the basis of bulk 
properties only. 

The paper is organised as follows: we introduce the 
model describing the magnetic properties of manganese 
oxide in Sec. II and study its properties in the clean limit 
in Sec. III. Next we discuss the mechanism which gen- 
erates the spin-glass phase in the doped samples and 
present the results of numerical simulations in support 
of our model in Sec. IV. Conclusions are presented at the 
end. 



II. MODEL 

The materials of interest for this study have a com- 
plicated lattice structure known as the hollanditc lattice. 
This structure consists of two octahedra joined at the 
edges to form the wall of 2 x 2 channels (see Fig. 1), 
and belongs to a family of crystals which differ from 
each other only by the lateral size of the channels, like 
ramsdellite (I x 2), romanechite (2 x 3) and todorokite 
(3 x 3)°' 24 ~ 2 '. In the case of manganese oxides each oc- 
tahedron is formed by a MnOg unit. 

The magnetism is due to the interaction of the mag- 
netic moments localised on the manganese ions, which 
interact with each other through oxygen-mediated su- 
perexchange 11,10 ' 28-30 . We consider the simplest possi- 
ble model compatible with the lattice structure and place 
classical Ising spins on Mn sites of the hollandite lattice. 
Therefore we get a lattice of spins which has 8 spins per 



unit cell as shown in Fig. 1. The Hamiltonian reads: 

n= J2 J l S i S 3+ X! J 2 S i S 3+ J2 J 3 S i S j, (!) 
<ij>l <ij>2 <ij>3 

where < ij >k denotes three different groups of nearest- 
neighbours of a spin. The partitioning of the neighbours 
and the corresponding couplings Ji, J2, J3 are detailed 
in Fig. 2. Such division and choice of the three different 
coupling constants are dictated by the the structure of 
the material: all these three classes of nearest neighbours 
are quite close to the atom: the distance to atoms of 
the 1 st group is ~2.86 A, while the atoms of the 2 nd 
and 3 rd groups are at distances of 2.91 A and 3.44 A, 
respectively. The choice of Ising spins is the simplest 
possible and, as we will show, it already reproduces the 
main experimentally observed features. 




Figure 2. (Colour online) The graphical representation of the 
Hamiltonian (1). Links of different colour represent the three 
couplings Ji (blue), Ji (black) and J3 (red). 



Despite the simplicity of the Hamiltonian (1) it admits 
many different phases depending on relative strengths 
and signs of the couplings Jk for k — 1,2,3. This is 
due to the complicated geometrical structure of the hol- 
landite lattice. If all the couplings are FM, the phase 
diagram only has paramagnetic and FM phases. Since 
the undoped compound is known to have an AFM or- 
dering experimentally—, we assume AFM couplings 
k = 1, 2, 3 throughout the rest of the paper. This choice 
corresponds to the maximum frustration of the interac- 
tions and reproduces experimentally observed features. 



III. CLEAN SYSTEM 

We start by studying the clean limit where the 
strengths of all of the couplings Jk k = 1, 2, 3 are constant 
and do not fluctuate in space (no disorder). Experiments 
indicate a transition to an AFM phase 22 at low enough 
temperatures. This behaviour is reproduced in our model 
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Figure 3. (a) (Colour online) C-type antiferromagnetic (C- 
AFM) ground state on the hollandite lattice at T — 
for 2J3 > J. There is an AFM alignment in the clock- 
wise/counterclockwise direction, (b) Staggered magnetisation 
m s as implied by the ground state of (a) as a function of tem- 
perature for the clean system (D=0), N = 2048 spins, volume 
V = 256Vb , where Vb is the volume of the unit cell. The cou- 
plings are equal to Ji = J2 = J = 1 and J3 = 0.6J. The 
inset: the staggered susceptibility Xs as a function of tem- 
perature. The sharp peak signalises the transition from the 
paramagnetic phase to Neel order. 



if we set the couplings to J\ = J2 = J and J/2 < J 3 < J 
where we have taken J = 1 for all the calculations. We 
assume J\ = J2 since the structural triangles formed by 
Mn atoms are almost equilateral 20 (see Fig. 2). With 
this choice of the interaction strengths there is a unique 
ground state with spins ordered in a C-typc AFM (C- 
AFM) state as shown in Fig. 3 (a). 

The C-type AFM order (C-AFM) is composed of FM 
chains which are coupled antiferromagnetically ,J (see 
Fig. 3). We now give a simple minimisation argument in 
favour of the C-AFM ordering for this particular choice 
of the coupling strengths. As can be seen from Fig. 2 
every Mn atom Sj interacts with 8 other Mn atoms 
Sik,k = 1,...,8. The coupling strengths for these 8 
atoms are as follows: two with J\ (sn,Si2), two with J2 
(si3,Si4) and four with J3 (si m ,m = 5,6,7,8). Therefore 
as long as 2 J3 > J, there is a unique state that minimises 
the total energy and it has s; m = — s,j = — s,-, j = 1,2. 
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This arrangement of spins does not minimise the energy 
for the links in the y direction sn = s,2 = Sj where spins 
are aligned ferromagnetically and the AFM links are un- 
satisfied. However this is exactly compensated by the 
interaction energy with the spins Si3 = Sa = — s,. The 
ground state average energy per spin is < Eq >= — 2J 3 . 
This ordered C-AFM configuration can be described as 
follows: if one numerates the 8 spins in the unit cell in 
the clockwise (or anticlockwise) direction s n ,n = 1, . . . , 8 
(see Fig. 3 (a)) then s n = (— 1)". This unit cell can be 
replicated on all the three directions implying that in the 
direction of the channels, along the y axis, the spins are 
aligned ferromagnetically (see Fig. 3 (a)). 

We have checked the validity of the above argument 
by performing classical Monte-Carlo simulations '" of a 
clean system with vV = 2048 Ising spins and a simula- 
tion cell of volume V = 256Vo where Vq is the volume 
of the unit cell composed by 8 spins. We assumed peri- 
odic boundary conditions and we have chosen J3 = 0.6 J, 
Ji = J2 = J- The standard Metropolis 32 rule with single 
spin-flip updates was used. We computed the staggered 
magnetisation which we defined as m s = M s /N with: 

N 

M s = £(-1)'*. (2) 

i=l 

where the site index i is the same as in the minimi- 
sation argument presented above and in Fig. 2 (ix_s)- 
The values of m s computed using the Monte-Carlo for 
various values of temperature are shown in Fig 3 (b). 
Clearly m s — for high temperatures and m s 1 
for temperatures close to zero indicating the AFM or- 
der. We also computed the staggered susceptibility Xs = 
dm s /dh s \ hs=0 = (< A/ 2 > - < M s > 2 )/k B T where 
ks is the Boltzmann constant and T the temperature. 
As one can see in the inset of Fig. 3 (b) the staggered 
susceptibility has a sharp peak at the same temperature 
where m s becomes non-zero. The critical temperature 
T c = 0.67 J/ks is estimated from the divergence of Xs 
(see inset of Fig. 3 (b)). 

For 2J3 < J the situation is completely different: the 
C-AFM phase is suppressed and the staggered magneti- 
sation is always zero: m a — as confirmed by Monte- 
Carlo simulations. Instead there is an extensive number 
of spin configurations that minimise the total energy as 
we show below. The ground states in this phase can be 
rationalised as follows: if 2J3 < J then the coupling to 
the nearest-neighbours in the 1 st and the 2 nd groups is 
more important than interaction with the 3 rd group of 
nearest-neighbours, unlike the case 2J3 > J. To better 
explain the spin configurations that minimise the total 
energy let's concentrate on one spin s, and assume Sj = 1 
(see Fig. 2). Since 2J3 < J, Sj interacts the strongest 
with the nearest-neighbour spins from the 1 st and the 
2 nd groups. Therefore we select sn = s;2 = — Sj = — f , 
and this creates an AFM order along the chain of Mn 
atoms in the y direction. Next we have to fix the spins 
Si3 and s,4, we notice that each of them interacts equally 



with two spins, one up and one down, in the chain con- 
taining Si (see Fig. 1 (c)), and therefore the interaction 
energy will be the same irrespective of the orientation 
of the spins 5^3 and s^. Then the energy is reduced 
by minimising the interaction between and Sn, that 
is setting s^ = — and creating a second AFM chain 
along the y direction that contains the spin. In this 
configuration the interaction of Sj with the sa and sn 
cancels out. Repeating the same reasoning for the other 
spins, we minimise the energy for the nearest neighbours 
of the I s * and the 2 nd groups in all the four Mn triangular 
chains (see Fig. 1 (b)). The interaction of Si of with spins 
of the 3 rd nearest-neighbour group also cancels out since 
Si5 = —Si6 an d = —Sis and the average total energy 
per spin is < Eq >= — J. To summarise this picture one 
can think that the ground state is composed by a collec- 
tion of one-dimensional Mn-Mn chains each having AFM 
order of spins. However the chains arc completely un- 
corrected between them and arranged in such a way to 
form the channels of the hollandite structure (see Fig. 4 
(a)). That is there is a perfect AFM order in y direction, 
but no order in the transverse direction. If we consider a 
single unit cell replicated in the y direction then we will 
have 2 8 possibles configurations with the same energy per 
spin < Eq > = — J. Since each chain is completely uncor- 
rected to the others, for an increasing number of cells in 
the x and z directions, A^j,the number of ground states 
is growing as 2 8Ar ™» , i. e. exponentially. 

The Monte-Carlo results confirm this picture: for 
2 J3 < J the staggered magnetisation is always zero and 
the staggered susceptibility does not display any sharp 
features down to T = 0. Such behaviour is typical 
for systems where the combination of lattice geometry 
and interactions suppresses the natural AFM order. At 
low temperatures, instead, they enter a phase known 
as collective paramagnet with power- law/exponential 
correlations depending on the properties of the lattice 
(bipartite/non-bipartite). We have computed the spin- 
spin correlation functions along different axes of the unit 
cell as shown in Fig. 4 (b) . The simulation was performed 
for the same system as presented in Fig. 3 (b) but with 
J3 = 0.3J. Unlike other geometrically frustrated sys- 
tems there is a clear AFM alignment of spins along the 
y direction and a (fast) exponential decay as we deviate 
from this direction (x, z 7 xy 7 xz,yz,xyz) and also with 
the 8 spins of the unit cell. This confirms the picture 
developed above, based on simple energy minimisation 
argument. 

For 2 J3 = J the system is even more degenerate: both 
families of the ground states introduced above have equal 
energies, allowing also for configurations that locally mix 
both solutions, making it difficult to put forward a min- 
imisation argument for this case. The system has an 
exponential number of ground states. It has no order at 
low temperatures and the AFM alignment along the y 
direction is suppressed. We believe that this lends sup- 
port to the view that the system might be in a classical 
spin liquid phase. 
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Figure 4. (a) (Colour online) Single unit cell replicated in the 
y direction. The spin configuration corresponds to one of the 
possible ground states with energy per spin < Eo >= — J. 
The spins are always aligned antiferromagnetically in the y 
direction for these ground states, while there is no order in the 
perpendicular direction, (b) The spin-spin correlations in the 
geometrically frustrated phase for the x,y, z,xy,xz,yz,xyz 
directions and inside the unit cell. Here the index is defined 
as the number of times you have to move in each direction 
to find Sj e.g. in the xy direction index=2 means that Sj 
is located at a distance of y (2L X ) 2 + (2L y ) 2 where L x and 
L y are the length of the unit cell in the x and y direction 
respectively. In the case of "Cell" the index runs over the 8 
atoms of the unit cell. 



To summarise: in the clean limit, with only AFM cou- 
plings there are three types of ground states. The one 
with C-AFM order (2J3 > J), a second one where there 
is a perfect anticorrelation along the one-dimensional 
Mn-Mn chains and no order is present between chains 
(2 J3 < J) and a third ground state consisting in a mix- 
ture of the previous two (2J3 = J). 



IV. DISORDERED SYSTEM 

We now turn to the disordered case and study what ef- 
fect the disorder has on the physical properties of the sys- 
tem described by the Hamiltonian (1). A spin-glass phase 
was experimentally observed in different compounds, for 
example KMngOi6 material 10 and in the case of Na + 
and K + doped a-MnOa 20 ' 21 ' 23 . Usually spin-glass be- 



haviour is associated with the presence of disorder and 
frustration in the system 33 . As discussed above, the clean 
system has geometrical frustration. We suggest the fol- 
lowing mechanism to explain the experimental results: as 
the dopant ions penetrate into the channels of the hol- 
landite structure 19 , they locally modify the compound 
and therefore also the magnetic interactions between the 
spins, i.e. they generate fluctuations in the AFM cou- 
plings Jfc, k = 1,2,3. In principle, this could happen 
through different mechanisms, such as doping-induccd 
strain or changes in the local electronic structure through 
charge donation. 

In the case of a-Mn02 it has been experimentally ob- 
served that lattice strains induced by doping are mi- 
nor 20 , while it seems more probable that the charge do- 
nated by the doping elements modifies the oxidation state 
of manganese, producing a local mixture of Mn +3 and 
Mn +9 ' 20 ' 34 . This local change in the electronic prop- 
erties induces fluctuations in the strength of magnetic 
spin-spin interaction, breaking the symmetry of the sys- 
tem (all Mn atoms are equivalent in the clean case). In 
our model we mimic these fluctuations by introducing 
continuous disorder in the couplings Jk. We would like 
to stress that we assume fluctuations of the couplings 
around their clean values so that the couplings remain 
AFM. The presence of quenched fluctuations in the cou- 
plings Jfe is crucial for the appearance of the spin-glass 
phase . Since the compounds investigated in the ex- 
periments have an AFM order at low temperatures in 
the clean limit and the appearance of the spin-glass re- 
quires doping to be higher than some critical value, this 
places us in the 2 J3 > J part of the phase diagram of our 
model. 

We propose the following mechanism for the appear- 
ance of the spin-glass in the doped compounds: random- 
ness in the AFM couplings Jk, k = 1..3 produces local 
regions of the lattice where the condition 2J3 > J is 
violated and that local patch has the C-AFM ordering 
locally disfavoured. The presence of such patches frus- 
trates the system, and generates the spin-glass phase. 
Such mechanism implies that the spin-glass can only ap- 
pear for sufficiently strong fluctuations in the strengths 
of the couplings: one needs sufficiently big number of the 
patches. Therefore we expect the critical strength of the 
disorder to be an increasing function of the distance to 
the point 2J3 = J. Since the amplitude of the fluctua- 
tions (disorder) depends on the doping, we need a suf- 
ficiently large doping to produce the spin-glass. This is 
exactly what is seen in experiments: Luo et a/. 20,21 have 
found that a critical amount of K + doping in a-K^MnC^ 
is needed to see the spin-glass phase. 

We have performed numerical simulations of the model 
with disorder to confirm our ideas. The AFM couplings 
Jij = Jji = Jk + A arc different for each link connecting 
a pair of atoms i,j = 1..N. The disorder A is drawn 
from the box distribution — D < A < D where D sets 
the scale for the strength of the quenched fluctuations. 
Glassy systems are notoriously difficult to simulate at 
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Figure 5. Staggered magnetisation m s as a function of tem- 
perature for the same system as in Fig. 3 except now a small 
amount of disorder D = 0.18 J drawn from the box distribu- 
tion is present in the AFM couplings. This disorder strength 
is not enough to destroy the C-AFM phase. Inset: staggered 
susceptibility \a as a function of temperature for the same 
system and the same disorder strength. 




Figure 6. (Colour online) Scaled Binder cumulant as defined 
in Eq. (5) for four system sizes V = 54 Vb, V = 96Vb, V = 
150Vb, V = 256Vb where Vb is the volume of the unit cell, 
the interactions Ji = J2 = J and J3 = 0.6 J where considered 
with a disorder of D = 0.4. Inset: Binder cumulant as defined 
in Eq. (3) as a function of temperature. 



low temperatures and we used parallel tempering to 
equilibrate the samples at low temperatures. We used up 
to N s = 200 samples to compute the disorder average. In 
every case we have tested that convergence to the average 
was achieved. 

We first show the results for small amounts of disorder 
D where the C-AFM phase is still present. In Fig. 5 we 
show the plot of the staggered magnetisation vs. tem- 
perature for D = 0.18J. The parameters of the sys- 
tem are the same as those of the clean case (D = 0, 
see Fig. 3):Ji = J 2 = J and J 3 = 0.6 J, N = 2048 
spins, V — 256Vo where Vq is the volume of the unit cell. 
We see that this amount of disorder is not enough to 
destroy the C-AFM phase but the critical temperature 
T c = 0.59J/fcs extracted from the position of the sharp 
peak of the staggered susceptibility (see Fig. 5) is reduced 
in comparison with the clean case where T c = 0.67 J/ks- 

Upon further increase of the strength of the disorder 
the C-AFM phase is completely suppressed at D = 0.4 
and the spin-glass phase appears instead at low tem- 
peratures. To characterize the paramagnetic to spin- 
glass transition we have studied the Binder cumulant 38 ' 39 
which is defined as: 



(3) 



where (• • ■ ) denote both thermal average for a given dis- 
order and the average over the disorder, and q is the 
overlap between two independent replicas of the system 
(sj, sf ) with the same realisation of disorder: 



1 N 



si 



(4) 



The Binder cumulant G is a dimcnsionlcss parameter 
that goes to zero for high temperatures and it is of order 
one in the spin-glass phase and has the finite size scaling 
form 3 " close to the transition: 



G = g \v 1/3 "(T - T c 



(5) 



where V is the volume of the sample and T c is the critical 
temperature. These properties make the Binder param- 
eter G a useful tool to study spin-glass phase transitions 
and it has been widely used for this purpose , as all 
curves of G generated for different system sizes will in- 
tersect at T c . As a complementary parameter we have 
also computed the overlap distribution P(q) that has the 
finite size scaling form : 

P(q) = V f)/ ^p{q) \qV^ 3u , V 1/3 »(T - T c ))l (6) 

where v and /3 are the critical exponents which are ob- 
tained by fitting the data for different system sizes to a 
single master curve. 

The parallel tempering simulations were performed for 
systems with AFM couplings J\ = J 2 = J, J3 = 0.6 J and 
disorder D = 0AJ. The system sizes are V = 54Vq,V = 
96Vb, V = 150Vb, V = 256Vb- In all the four cases the 
C-AFM phase was suppressed in favour of a spin-glass 
phase. The inset of Fig. 6 shows the Binder cumulant 
as a function of temperature for the four system sizes. 
There is a clear crossing of the Binder cumulant curves 
indicating the spin-glass transition near T c = 1.19J/&B 
for D = 0.4J. 

Using this estimate for T c we collapsed all the data onto 
a single master curve according to Eq. (5) as shown in 
Fig. 6 whit the critical exponent v = 0.86. We have also 
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Figure 7. (Colour online) Scaling for the probability distri- 
bution P(q) at T = 1.20J/fc_g near the critical temperature 
T c = 1.19J/fc_e using Eq. (6) for the same system of Fig. 6. 



studied the overlap distribution P(q) for the four system 
sizes at the temperature T = 1.20J/fcs, close to the crit- 
ical temperature T c = 1.19J/fcs (see Fig. 7). Rescaling 
the P(q) curves for different system sizes according to 
Eq. (6) a good collapse on a master curve is obtained for 
the ratio of the critical exponents j3/v = 0.42 from witch 
we get (3 = 0.36. 

The values j3 = 0.36 and v = 0.86 are very different 
from the values of the critical exponents for the spin-glass 
transition in d = 3: /3 = 0.77(5), u = 2.45(15) (see Ref. 44 
and references therein). These exponents are known to 
be hard to evaluate precisely due to large finite-size ef- 
fects 44 . As there are no other factors that could explain 
the discrepancy we attribute it to the smallness of the 
system sizes studied. Understanding better this discrep- 
ancy and studying in detail the dependence of the critical 
exponents on the system size requires to simulate much 
larger system sizes which is beyond the scope of this pa- 
per. 



V. CONCLUSIONS 

We have introduced a simple classical Ising spin model 
with nearest-neighbour antiferromagnetic (AFM) inter- 
actions and continuous disorder. This model explains 
recent experimental results on the magnetic properties 
of some manganese oxides with hollandite crystal struc- 



ture. The studied system is geometrically frustrated due 
to the existence of structural triangles (see Fig. 1 (c)) 
in the hollandite lattice. The parallel tempering Monte- 
Carlo simulations performed for this system with clean 
and disordered AFM coupling strengths and different sys- 
tem sizes, reveal that, despite the simplicity of the Hamil- 
tonian Eq. (f )) (only nearest neighbour interactions), it 
admits different phases depending on relative strengths 
of the AFM couplings due to the complicated geometry 
of the hollandite lattice. In fact there are three types 
of ground states in the clean case: a non-degenerated 
one with a C-type AFM (C-AFM) order, composed of 
FM chains coupled with an AFM order between them 
(see Fig. 3). The other two types of ground states are 
exponentially degenerated and disordered ground states. 
Depending on the strengths of the interactions a small 
amount of disorder in the couplings does not suppress 
the C-AFM order, but merely reduces the critical tem- 
perature at which the paramagnct to C-AFM transition 
occurs (see Fig. 5). When the disorder becomes greater 
than some critical value, the geometry forces different 
ground states to compete and there is a direct transition 
from the C-AFM state to a spin-glass by tuning the dis- 
order (at fixed temperature) , while preserving AFM cou- 
plings. Such a transition was confirmed in numerically 
by looking at the Binder cumulant (see Fig. 6). To the 
best of our knowledge this is the only existing model that 
presents this kind of transition with nearest-neighbours 
AF interactions. 

This model represents a good starting point to under- 
stand the magnetic phases appearing in hollandite oxides. 
It reproduces the main features of the experimental re- 
sults both in clean and doped cases. This model shows 
that bulk effects (frustration and doping) are sufficient to 
obtain the spin-glass in compounds with hollandite lat- 
tice. It would be interesting to check our predictions ex- 
perimentally and to see how adequate the proposed mech- 
anism for generating the spin-glass order is. It would also 
be useful to see the convergence of the critical exponents 
to their standard d = 3 values, although this requires 
considerable amount of computing time. Finally the pre- 
dictive power of our model could be enhanced by using 
couplings determined by ab-initio calculations. 
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